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A measurement of the ratio R^a ~ iXc -^ J/i' + 1)1 J !'']> in pC, pTi and pW interactions at 920 
GeV/c (y^ = 41.6 GeV) in the Feynman-x range —0.35 < Xp < 0.15 is presented. Both /x "*"//" 
and e+e^ ,J/ip decay channels are observed with an overall statistics of about 15000 Xc events, which 
is by far the largest available sample in pA collisions. The result is R-,^^ = 0.188 ± 0.013stlgQ22 
averaged over the different materials, when no J/i/> and Xc polarisations are considered. The Xci to 
Xc2 production ratio R12 = Rxail Rxa2 ^^ nieasured to be 1.02 ±0.40, leading to a cross section ratio 
^^'^"^j = 0.57 ±0.23. The dependence of R^^ on the Feynman-x of the J/ip, xj, and its transverse 

momentum, Py , is studied, as well as its dependence on the atomic number, A, of the target. 
For the first time, an extensive study of possible biases on R-^^ and R\2 due to the dependence 
of acceptance on the polarization states of J/ip and Xc is performed. By varying the polarisation 
parameter, Xohs, of all produced J/xp's by two sigma around the value measured by HERA-B, and 
considering the maximum variation due to the possible Xci and Xc2 polarisations, it is shown that 
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13.85.-t Hadron-induced high- and super-high-energy interactions 

24.85.+P Quarks, gluons, and QCD in nuclei and nuclear processes 
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I. INTRODUCTION 

Since the discovery of charmonium more than thirty 
years ago, its production in hadronic coUisions has at- 
tracted considerable theoretical and experimental inter- 
est for a variety of reasons. In particular the ques- 
tion of the production mechanism, which requires an 
understanding of the hadronisation process in the non- 
perturbative regime and, in addition, the influence of 
nuclear matter, are of particular importance since the 
suppression of J /ip production has been considered as a 
possible indicator of the quark-giuon plasma []| . 

The theoretical treatment of quarkonium production 
is usually broken into two steps: the creation of a 
heavy quark pair in interactions of the colliding par- 
tons, calculable by means of perturbative QCD, and 
the transition to a bound state, involving poorly under- 
stood non-perturbative processes and even more prob- 
lematic nuclear effects. A variety of approaches have 
been developed to describe quarkonium production such 
as the Color Evaporation Model (CEM) [2], the Color 
Singlet Model (CSM) 3, and non-relativistic QCD 
(NRQCD) (^j. A measurement of the fraction of J/tp 
coming from the decay of other charmonium states (feed- 
down) provides useful tests of the model predictions. 
While a rather rich sample of data on J/^ production ex- 
ists, the available data on the production rates, or even 
the experimentally simpler fractional production rates, 
of the other charmonium states suffer from imprecision. 
Moreover very little experimental information is available 
on the possible polarisation of the produced charmonium 
states. 

In this paper we report on the production of the char- 
monium states Xci and Xc2 in collisions of a 920 GcV 
proton beam with nuclear targets. The Xc mesons are 
identified via their radiative decay into J /ij) mesons which 
in turn are decaying into lepton pairs. The production 
and decay chain is: 



pA^Xc + X; 



Xc 



-fJ/^^-fl+l- a = e, Ai). (1) 



To minimise systematic uncertainties the Xc rates are 
normalised to the total production rate of J/ip. We define 
RxoJ the fraction of J/ip originating from radiative Xc 
decays: 



Ry 



E^=l (^{Xc^)Br{Xc^ "> J/^l) 



(2) 



where Br{xci ~> J/ipl) a-re the branching ratios for the 
different Xci -^ J I'Pl decays, cr{Xci) are their production 
cross sections per nucleon and a{J /ip) is the total J/i/i 



production cross section per nucleon. In Tab. [J the main 
properties of the three Xc states (xcOi Xci and Xc2) are 
reported. Due to the negligible Xco -^ J l'4'l branching 
ratio, we limit our study to Xci and Xc2 production. 

State Mass (MeV/e^) Width (A/eV/e^) BR(^ J/'ip-y) 



Xco 3414.76 ± 0.35 
Xci 3510.66 ± 0.07 
Xc2 3556.20 ± 0.09 



10.4 ±0.7 
0.89 ±0.05 
2.06 ±0.12 



(1.30 ±0.11)% 
(35.6 ±1.9)% 
(20.2 ±1.0)% 



TABLE I: Properties of the three Xc states 0/. 



In Scction|TT]an overview of the physics motivations for 
our measurement is given together with a survey of the 
existing experimental results. In Section IIIII the experi- 
ment and the data taking conditions are described and in 
Section HVl the principle of the measurement is explained. 
The simulation used for evaluating detection efficiencies 
and the measurements in the muon and electron decay 
channels are described in Sections |Vl IVI Al and IVIBi re- 
spectively. In Section IVIIi the effect of possible J/^ and 
Xc polarisations (not directly measured in this analysis) 
on R^^ and R\2 = Rxai/Rxa2 ^^ discussed. A discussion 
of systematic uncertainties and the final result are given 
in Sections IVIIII and IIXI respectively. 



II. Xc PRODUCTION 

A. QCD Models of Charmonium Production 

In the Color Evaporation Model (CEM), charmonimn 
production is described as the creation of cc pairs with 
an invariant mass below the DD threshold. Their hadro- 
nisation is mediated by the emission of soft gluons which 
do not significantly alter the kinematics of the cc sys- 
tem. The production rates of the various charmonium 
states are predicted to be proportional to each other and 
independent of the projectile, target and energy. Most 
experiments measuring R^^ in proton and pion induced 
interactions [6|-[2l| provide compatible results with the 



predicted value R^ 



0.4 



The assumption of the 



universality of charmonium hadronisation implies that 
the value of R^^ should be independent of the kine- 
matic variables xp and px of the produced charmonium 
state [22| {xp is the Feynman variable in the nucleon- 
nucleon center-of-mass system; pT is the transverse mo- 
mentum relative to the incoming beam). 

In the Color Singlet Model (CSM), the quark pair is 
created in a hard scattering reaction as a colour sin- 



glet (CS) with the same quantum numbers as the final 
quarkonium. Since two gluons can form a colourless C- 
even state such as the Xc states, but at least three gluons 
are needed to form a colourless C-odd state such as the 
V' states, TJ) production in this model is suppressed by an 
additional factor Us- As a result, the J/-0 production 
rate should be dominated by feed-down from radiative 
Xc decays and R^^ is predicted to be close to 1. Most of 
the proton induced Xc measurements are in disagreement 
with this assumption |6l|-|14|. 

In response to the disagreement between the CSM 
and measurements in most charmonium production 
features [23|, a more generahsed perturbative QCD 
approach for charmonium production, 'non-relativistic 
QCD' (NRQCD), was developed which includes not only 
cc pairs produced as colour singlets but also as colour 
octets (CO). The CO states subsequently evolve into the 
observed charmonium by soft gluon emission. At the 
HERA-B beam energy of 920 GeV, the dominant pro- 
duction process is gg fusion which contributes both to 
CO and CS states. Therefore Xc production dominates 
the CS part of J/-0 production while direct J/-0 and J/-0 
from ^jj' decay are produced via CO states. The pre- 



dicted ratio, R^ 



0.3 



is in agreement with most 



of the exis ting measurements in proton induced inter- 
actions [y|-[ij]- NRQCD predicts only small differences 
in the differential cross sections of the different charmo- 
nium states as a function of xp, mostly at large values of 
xp- More visible differences can arise when considering 
nuclear-matter effects (A-dependence) due to the differ- 
ing absorption probabilities of the various charmonium 
and pre-charmonium states in nuclei [22j . 



B. Interactions with Nucleons and A-dependence 

The CEM model and NRQCD differ in their predic- 
tions of the suppression of the charmonium production 
rate per nucleon in interactions with heavy nuclei com- 
pared to interactions with single proton targets. Sup- 
pression can occur in interactions of the generated cc 
quarks with nuclear matter which could lead to an xp 
dependence: for xp > 0, the formation length of the final 
charmonium state exceeds the size of the nucleus, while 
for xp < 0, an increasingly larger fraction is formed al- 
ready inside the nucleus. In the context of the CEM, only 
one proto-charmonium state exists and thus for xp > 0, 
the differences in suppression between J/f/', V'' ^nd the 
Xc states should be small. For xp < 0, the Xc and ip' 
states should be more suppressed than the J/iJj due to 
their larger interaction cross sections. In the context 
of NRQCD, substantial differences in the suppression 
of the various charmonium states are expected even for 
Xp > since the wave function of the CO states extends 
over a much larger distance and the resulting interaction 
cross section is considerably larger than that of the CS 
states [22 1 . 



C. Ry.^ and the Quark-Gluon Plasma 

The so-called "anomalous" suppression of J/ip has 
been proposed as a possible indicator of the forma- 
tion of a quark-gluon plasma [l|] and such suppres- 
sion has subsequently been reported by several experi- 
ments |2J| [2^ [201 . Nevertheless the conclusion that the 
reported suppression is indeed anomalous is contingent 
on the full understanding of normal suppression mech- 
anisms, i. e. those existing in the absence of a quark- 
gluon plasma as is expected to be the case in proton- 
nucleus reactions. In this respect the measurement of 
the fraction of J/tp arising from feed-down decays {xc 
and tp') is important since the anomalous suppression 
is expected to be sensitive to the mass and binding en- 
ergy of the different charmonium states. Directly pro- 
duced J /lb survive in the quark-gluon plasma up to about 
1.5 Tc |27| . Tc being the critical temperature, while Xc 
and ip' states dissociate just above Tc- Thus several drops 
in the distribution of charmonium survival probability 
as a function of the temperature are expected, with the 
size of the drops dependant on the fractions R^^ and R^i 

{R^' = '^ ct(j/^ )• Experimentally, only the first 

drop has been reported [241 [2a| [2a] , and is interpreted as 
indicating the disassociation of Xc and tp' . Several mod- 
els attempt to describe the totality of experimental data 



on J/^ suppression. They generally assume R^ 



0.3 



and R^' ~ 0.1. Nevertheless all the proposed models fail 
to simultaneously describe all the existing data, as they 
all overestimate the suppression unless other effects, such 
as J/V^'egeneration, are assumed to describe the RHIC 
data [23. From the value of R^^^ shown in Sect. IIXI and 



the result from '28|, i?^- « 7%, R^^ 
is lower than generally assumed. 



+ R> 



■Ip' 



0.27, which 



D. Previous Measurements 

The production of Xc has been measured both in 
proton- and pion-induced reactions on various nuclear 
targets and in pp and pp interactions [a|-[21|. Tab. [llj 
lists all the published measurements of Xc production in 
hadronic interactions and reports their most relevant fea- 
tures. From this table some observations can be made: 

• all fixed target measurements are based on at most 
a few hundreds Xc] 

• all experiments observe only one of the two J /tp 
decay channels (e~''e" or /i+/x^); 

• the photon efficiency never exceeds 30%; 

• most measurements are performed in the positive 
Xp range. 

Tab, mil shows the measured values of R^^ and/or of the 
Xc cross sections separately for proton and pion induced 
reactions. The values shown in Figs. [T] and [2] have been 



updated using the current PDG values \3\ for the Xc and 
J/ifj decay branching ratios, and the J/ip cross sections 
obtained from |29j . 

The available data scatter strongly, well beyond their 
respective uncertainties, and no energy dependence is 
discernible. The proton data seem to favour a value 
R^^ ^ 0.3, supporting the prediction of NRQCD, but 
the quality of the available data does not allow a firm 
conclusion. 



of the photon from Xci (Xc2) decay for events with iden- 
tified J/-0's. The measurement method consists of eval- 
uating Nj/jf, by analysis of the di-lepton invariant mass 
spectra and N^^-^ , N^^^ by analysis of the J/ip-J invariant 
mass spectra, for events with selected J/ip candidates. 
The efficiency terms in Eq. [3] are extracted from the MC 
simulation. 

The production ratio of the two states can be deter- 
mined using: 



III. THE EXPERIMENT AND THE DATA 

SAMPLE 

The HERA-B detector [30| was a forward magnetic 
spectrometer used to study the interactions of the 920 
GeV proton beam (^i = 41.6 GeV) of the HERA ac- 
celerator on a variety of nuclear targets. The detector 
components relevant for this analysis are the wire tar- 
get system [31] which could be dynamically positioned 
in the halo of the proton beam, the Silicon Vertex De- 
tector (VDS) [13], the dipole magnet of 2.13 Tm, the 
drift-tube Tracking System (OTR) [13], the Ring Imag- 
ing Cherenkov Counter (RICH) 533], the sampling Elec- 
tromagnetic Calorimeter (ECAL) [35] and the Muon De- 
tector (MUON) 1361]. 

The data sample of about 160 million events used for 
this analysis was acquired at an interaction rate of about 
5 MHz with a dedicated di-lepton trigger [33] in order 
to select both J/ip — > e+e~ and J/tp — > /i"'"/i~ final 
states. In total about 300,000 J/i/' were reconstructed, 
distributed almost equally in the two decay channels. 
Nine different wire configurations were used, both in 
single and double wire runs. The wire materials used 
were carbon (C, «64% of the full statistics), tungsten 
(W, «31%), and titanium (Ti, «5%). Continuous on- 
line monitoring ensured stable running conditions, and 
further offline data quality checks were applied to select 
only runs with properly functioning detector and trigger 
components. 



IV. EXPERIMENTAL METHOD 



R^^ is defined in Eq. [21 The quantity to be measured 



is: 



N^ 



Afv 



Ry 












N 



iVv 



J/i> 






1) 



Wv 



-J/± 
_Xc2 



1) 



(3) 
where iVj/^ is the total number of observed J ftp's, N^^-^ 

(-^Xc2) is the number of counted Xci's (Xc2's), Ej'L is 
the direct J/tp total detection efficiency, including trig- 
ger losses, reconstruction and cut selection, 6^",^ (cj/l) 
is the total detection efficiency for J/tp coming from Xci 
(Xc2) decay and e^"'^ i^^"^) is the identification efhciency 



Ri2 — 



R 



Xal 



AT F^"^ ■ F^"^ 



R 



-Xcl ^Xcl 



Xc2 



^Xc2 £j/'^ • £7 



(4) 



(where i?xci + ^Xc2 ~ -^Xc) and the production cross 
section ratio can be evaluated using: 



g(Xcl) _ p ^''(Xc2 

—7 T - -K12 

C^(Xc2) 



J/ijj) 



Brixci -^ J/iPl) 



(5) 



In order to perforin an internally consistent analysis, 
the same procedure and cuts are applied to both the e+e~ 
and /i^/i^ channels except for the lepton particle identi- 
fication (PID) requirements. 



A. J/tp Selection 

Leptons from J/tp decay are selected from the trig- 
gered tracks, re-fitted using offline alignment constants 
and taking into account multiple scattering when extrap- 
olating to the target. A x^ probability of the track fit 
> 0.3% is required. Additional PID cuts are applied de- 
pending on the lepton channel. 

In the muon channel, a muon likelihood is constructed 
from information in the MUON detector and is required 
to be greater than 5% and a kaon likelihood is con- 
structed from RICH information and required to be less 
than 99%. 

In the electron channel a more complex set of PID cuts 
is needed. First the calorimeter is searched for a cluster 
consistent with having been caused by a Bremsstrahlung 
photon emitted in front of the magnet [3_5]. Since the 
presence of such a cluster very effectively identifies elec- 
trons, the cut values used for the remaining two particle 
identification criteria can be substantially relaxed when 
such a Bremsstrahlung cluster is found. The additional 
two criteria are a more restrictive matching requirement 
between the OTR track of the electron candidate and 
its corresponding ECAL cluster, and a requirement that 
the track momentum be consistent with the deposited 
calorimeter energy. 

Once opposite sign lepton candidates {fi'^fj,~ or e+e^) 
are selected, their common vertex is fitted and the x^ 
probability of the fit is required to be greater than 1%. 
In a few percent of the events, more than one di-lepton 
combination pass all cuts, in which case only the one with 
the lowest product of track-fit x^ is retained. Finally, 
the invariant mass of the di-lepton pair is calculated and 
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TABLE II: Previous R^^ measurements in hadronic collisions. Symbols: 7 detection (d— direct, c— ^-conversion). Xci-Xc2 
separation (y—yes, n=no, f=with 2-states fit). 



Exp. 



Measured values 



Updated values 



RXc 



g(Xcl) 
ct(Xc2) 



[6] 0.43 ±0.21 

m 



o-CXci) 
(nb/n) 



C^(X=2) 

(nb/n) 



RXc 



ct(Xc2) 



'^(Xcl) 

(nb/n) 



'^(Xc2) 

(nb/n) 



[8] 
[9] 
[10] 
[12] 
[14J 



l^s+oio 

U-J^O-O.lS 

0.47(8) 

0.47(23) 

0.30(4) 

0.77(30)(15) 

0.32(6)(4) 



0.24(28) 64(81) 268(136) 

0.08(25)(15) 31(62)(3) 364(124)(36) 
0.53(20)(7) 526(138)(64) 996(286)(134) 



0.43 ±0.21 

it^+oio 
'J-J^O-o.is 

0.47(8) 

0.47(23) 

0.30(4) 

0.76(29)(16) 

0.32(6) (4) 



0.24(28) 39(49) 162(81) 

0.09(29)(17) 24(48)(2) 244(83)(16) 
0.61(24)(4) 488(128)(56) 805(231)(92) 



[11] 0.297(17)(57) 1.04(29)(12) 



0.297(17)(57) 1.19(33)(14) 



[15] 
[16] 
[17] 
[18] 
[19] 
[20] 
[20] 
[21] 



0.70(28) 
0.30(5) 
0.44(16) 
0.37(9) 
0.31(10) 
0.40(4) 
0.37(3) 
0.443(41) (35) 



0.68(28) 

l(fix) 
1.12(42) 
0.96(64) 

0.70(15) 
0.57(16) 



65(18) 
28(10) 

130(56) 

131(17) 
464(87) 



96(29) 
28(10) 

134(64) 

189(31) 
815(168) 



0.70(28) 
0.30(5) 
0.44(16) 
0.37(9) 
0.31(10) 
0.40(4) 
0.37(3) 
0.443(41) (35) 



0.79(28) 

l(fix) 
1.11(41) 
0.98(74) 

0.80(16) 
0.65(18) 



58(13) 
22(8) 

102(43) 

101(13) 
356(66) 



74(19) 
22(8) 

104(49) 

126(19) 
544(107) 



TABLE III: Ry 



•-xc; ^()> ''"(Xci) o.'nd o{xc'2) results in hadronic collisions. Statistical and systematic uncertainties are shown 
in brackets (less significant digits). See text for an explanation of the updated values. 



required to be within 2a of the nominal J/^p mass, with 
cr = 36 MeV/c^ in the muon channel and 64 MeV/c^ in 
the electron channel. 



B. Mass Difference Plot 



The next step after J/ip selection is the identification 
of suitable photon candidates. A photon is defined as a 
reconstructed ECAL cluster [33| with at least three con- 
tiguous hit cells. The cluster energy, E"^ , is required to be 
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FIG. 1: Summary of experim,ental results on R^^ (a), "Jf^'^^l (b), (t(xci) (c) and a{xc2) (d) in pN interactions. HERA-B 
results: open star (old), full star (new). 
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FIG. 2: Summary of experimental results on R-^^ (a), '^f"'') (b), cr(xci) (c) and a{xc2) (d) in tvN interactions 



ct(Xc2) 



at least 0.3 GeV and the cluster transverse energy, E'^, 
is required to be at least 0.2 GeV, for an optimal clus- 
ter reconstruction. Furthermore, the ECAL cell with the 
highest energy deposit of the cluster is required to contain 
at least 80% of the total cluster energy in order to provide 
some discrimination against showering hadrons. Clusters 
which match reconstructed tracks are excluded unless 
the matching track is formed only from hits behind the 
magnet and point to the selected di-lepton vertex. Such 
tracks are mainly from conversions of event-related pho- 
tons behind the magnet. Finally, because of high back- 
ground near the proton beam pipe, clusters in an ellip- 
tic region around the pipe W ^tiust I '^ + v'tiust < 22 cm, 
where Xdust and Vciust are the horizontal and vertical 
positions of the cluster with respect to the beam) are 
excluded. 

Since a photon from a Xc decay cannot be distinguished 
from the others in the event (on average ~ 20), the com- 
binatorial background to the Xc signal is very large, as 
will be shown in Sect. IIVCI 

To largely eliminate the uncertainty due to di-lepton 
mass resolution, the analysis is performed using the mass 
difference, AM = M{J/tp-f) - M{J/ip). The domi- 
nant contribution to the mass difference resolution is 
the intrinsic photon energy resolution determined by the 
ECAL. 



C. Background Description 

The analysis crucially depends on the background 
shape being correctly described. We distinguish between 
"physical" backgrounds (due to the decay of heavier 
states which include a J/V' and one or more photons in 
their decay products) and "combinatorial" background 
(due to photons from the event combined with di-leptons 
which share no parent resonance). The combinatorial 
background by far dominates. The only significant phys- 
ical background comes from tp{2S) -^ J/?/'7r°7r° which 
contributes at the level of « 15% of the Xc rate but with a 
rather flat distribution in the AM spectrum. The shape 
of this background is estimated from Monte Carlo and 
subtracted after proper normalisation. 

A "Mixed Event" (ME) procedure is adopted for mod- 
elling the combinatorial background: a Jf-ip candidate 
from one event ( "event- A" ) is mixed with the photons 
of several (w 20) other selected events (which we all call 
"event-B"). Event-B is required to have the same neutral 
cluster multiplicity as event-A to ensure similar photon 
energy spectra. Furthermore, the angular difference be- 
tween the vector sums of transverse momenta of all pho- 
tons in event-A and event-B is required to be no more 
than 27r/20 to ensure the events to be kinematically sim- 
ilar and thus to have similar acceptance. 



Extensive tests, both with Monte Carlo and the data 
itself, were performed to verify the ME procedure. For 
example, using the data, the combination of photons with 
l~^l~ pairs in the J/i/j side bands (defined as the di-lepton 
mass intervals outside Scr of the nominal J/ip mass, see 
Sect. lIVA| in the SE ("same event") and with l+l^ pairs 
inside the J/f/; mass window in the ME spectra show no 
reflection of the Xc peak and the SE over ME ratio for 
these events is found to be flat. The normalisation of 
the ME spectrum is incorporated into the fit of the AM 
spectrum as a free parameter (see Sect. lVIBJ) . 



V. THE MONTE CARLO SIMULATION 

A. Event generator and detector simulation 

In the HERA-B Monte Carlo, the basic process pN -^ 
QQX is simulated, first, by generating the heavy quarks 
(QQ), including hadronisation, with PYTHIA 5.7 [39[; 
secondly, the energy of the remaining part of the process 
(X) is given as an input to FRITIOF 40], which is used 
to simulate the interactions inside the nucleus. PYTHIA 
describes by default the charmonium production based 
on the Color Singlet Model. Further colour singlet and 
colour octet processes were therefore added, according to 
the NRQCD approach [1J|. Differing kinematic distri- 
butions for directly produced J/tp and J/tjj from feed- 
down decays generated according to this model result in 
slightly different acceptances: ^ 78.3% for direct J/tp 
and ~ 77.6% for J/tJj from both Xci and Xc2, with no 
significant difference between the two Xc states. 

In the simulation, both direct J/ip and Xc states 
are generated with no polarisation and all results are 
given under this assumption. The effects of J/ip and 
Xc polarisation are discussed and treated separately (see 
Sect.lVHI. 

The detector response is simulated using GEANT 
3.21 [4l| and includes individual detector channel reso- 
lutions, noise, efficiencies and calibration precision. The 
second level trigger algorithm is applied to the simulated 
detector hits and the first level trigger efficiency is taken 
from an efficiency map obtained from the data itself. The 
generated Monte Carlo is reconstructed with the same 
package used for reconstructing the data and the same 
analysis cuts are applied to the MC and the data. 

In order to check the MC material description, 
which influences the photon efficiency determination, 
three different studies were performed by using the 
Bremsstrahlung tag [3a|, the tt" signal (where the de- 
cay photons are seen as neutral clusters or as converted 
photons), and the converted photons (see Sect. IV Bl 
andEIIII- 

The predicted resolution of the Xci a-nd Xc2 states is 
found to be ~ 0.032 GeV/c^, in agreement with real data 
(see Sect. EH]). 



B. J/tp and photon efficiency 

According to Eq. [31 the ratios of efficiencies for J/tp 
coming from the decay of the Xc states to that of di- 
rectly produced J/tp are needed. These ratios are esti- 
mated from MC and the values obtained are reported in 
Table IIVI As can be seen from the table, these ratios 
are independent of target, decay channel and Xc state, 
within the errors. 

The photon detection efficiencies are also evaluated 
with the MC although an additional correction factor 
derived from the data was found to be needed, as will 
be discussed below. For the efficiency evaluation, the 
same analysis as for the data is performed, but the pho- 
ton from the Xc decay is selected using MC generation 
information and checked for acceptance after all cuts are 
applied. The alternative of extracting the number of XcS 
from the MC using the ME background subtraction ap- 
plied to the data, and thus inferring the photon efficiency 
without recourse to the MC generation information, was 
found to give a compatible efficiency, but with lower pre- 
cision. 

The MC estimate for photon detection efficiency was 
checked by comparing the efficiency derived from MC to 
that, obtained from data, for the detection by the ECAL 
of reconstructed electrons or positrons from photon con- 
versions before the magnet. Since the average di-lepton 
triggered data run contains several thousands of such re- 
constructible conversions, the method affords a detailed 
check of the stability of photon detection efficiency over 
the run as well as a check of the MC. 
The tracks from the converted photons are required to 
share a common VDS track segment and to have hits 
in the OTR chamber immediately before the ECAL (to 
discriminate against electrons which start to shower be- 
fore the ECAL). When using the positron from such a 
pair as a probe, the electron ("tag") is also required to 
have an associated ECAL cluster with a deposited en- 
ergy compatible with the electron track momentum (and 
vice versa). For selected electron and positron probes, 
the ECAL is searched for a geometrically matching clus- 
ter and the ratio of the deposited ECAL energy to the 
track momentum ("E/p ratio") is entered into a his- 
togram. Signal to background ratios of the order of 15 
are achieved. The E/p ratio histogram of the probe as 
well as the corresponding E/p histogram of the tag are 
fitted to gaussians to describe the signal and third or- 
der polynomials for the background description. The fit 
describes the data well with x^ values typically equal to 
or less than the number of degrees of freedom. The effi- 
ciency is extracted from the fit parameters. The ratio of 
MC efficiencies estimates to the efficiency derived by this 
method is found to be 1.144± 0.034, with the quoted un- 
certainty dominated by run to run variations. Roughly 
half the difference between efficiency estimates from MC 
and data can be attributed to a higher ECAL cluster 
multiplicity in the data compared to the MC - when a 
cluster caused by a photon from a Xc overlaps with an- 



other cluster the photon's measured energy becomes too 
large and the mass estimate incorrect. The remaining 
(« 7%) discrepancy is not understood but is likely due 
to cases where the energy deposited by the photon (elec- 
tron) is considerably less than would be expected from 
gaussian statistics. 

In Tab. |V|the values of efficiency and the width of Xci 
and Xc2 are reported for the two lepton channels and for 
the different target materials. 



Mat. 



-3r 



M> 












c 


0.972(7) 


-Ji^ ^.iii, "j/i, 
0.965(5) 0.970(12) 0.950(7) 


w 


0.957(8) 


0.974(6) 0.985(14) 0.955(9) 


Ti 


1.008(26) 


0.957(17) 



TABLE IV: The ratio of efficiencies for detection of directly 
produced J j^l) 's to that of J/tp 's coming from xa and Xc2 
decay. The efficiencies for fi'^ fi~ ande'^e' channels and for 
each target material are given separately. 



Mat. 


£^- (%) 


a^,, (MeV/c^) 


£^- (%) 


^x.2 (MeV/c^) 


C 


40.5 ± 0.4 


30.2 ± 0.4 


41.2 ±0.2 


33.0 ±0.2 


W 


37.1 ±0.5 


32.0 ±0.7 


38.2 ±0.3 


34.3 ±0.5 


Ti 


41.3 ±1.3 


31.6 ±0.8 


41.4 ±0.8 


30.2 ±0.3 


C 


39.6 ±0.6 


32.1 ±0.4 


40.4 ±0.3 


33.0 ±0.2 


W 


38.6 ±1.0 


33.7 ±0.8 


38.3 ± 0.6 


35.4 ±0.4 



TABLE V: Photon detection efficiencies and the expected 
widths of AM peaks for Xci and Xc2 for the muon (first part) 
and the electron (second part) channels. 



VI. EVENT COUNTING 

A. J/tp counting 

a. The muon channel: The /i^/x~ invariant mass 
spectra for C, Ti, and W samples as well as the summed 
spectrum are shown in Fig. [3] along with a fitted curve. 
The fit includes the J/4> and tp' peaks, each described by 
a superposition of three gaussians with a common mean 
plus a radiative tail to describe the photon emission pro- 
cess ip -^ iX^^~^ [2^, and an exponential to describe 
the background. The numbers of J/-0 within the mass 
window used for Xc selection are reported in Tab. IVll 

h. The electron channel: The e"'"e~ invariant mass 
spectra for the different materials and the full sample are 
shown in Fig. [H The fit used for the signals {J /ip and 
V^') includes a Gaussian for the right part of the peaks 
and, for the left part, a Breit-Wigner to take into account 
the Bremsstrahlung tail, while a gaussian (exponential) 



describes the background in the low (high) mass region 
with the requirement of continuity of the functions and of 
the first derivatives. The numbers of J /tj} within the mass 
window used for Xc selection are reported in Tab. I VII 



B. Xc counting 

The AM spectra are shown in Fig. [5] for the muon 
channel, and Fig.lHlfor the electron channel. In the upper 
parts of these figures, the SE data arc indicated by points. 
Fits to the ME and SE samples are shown as solid lines. 
The two curves are not distinguishable except in the AM 
region between 0.3 and 0.6 GeV/c^ where the ME curve 
is below the SE curve. The fit to the SE spectrum uses 
the ME parameterization to describe the background and 
Gaussian distributions for the signal, as described below. 
In order to evaluate the quality of the background de- 
scription, the background subtracted spectra are shown 
below the fitted AM spectra for visual representation 
only. A clear Xc signal is visible both in the carbon and 
tungsten samples, while in the titanium sample, the sig- 
nificance of the Xc signal is at the level of ~ 3.5tT only. 

The detector resolution for the two Xc states is compa- 
rable with their mass difference, resulting in a single Xc 
peak in the AM spectrum. It is nevertheless possible to 
count separately the number of Xci and Xc2 (and there- 
fore measure R12) by using a fit with two gaussians with 
some of the parameters fixed, namely: 

1. AM^^, = 0.4137 GeV/c2 [|; 

2. CT^^j fixed according to the MC prediction from 
Tab El 

3. AM(xc2) - AM(xci) = 0.0455 GeV/c^ [|; 

4. ^ = 1.05 as predicted by MC. 



^Xcl 



^V 



TV 



^^sJ- and the 



The free parameters are N^^ — jv^^^+^^j, j^ 
ME normalisation parameter. In Tab. IVIII the values of 
the fitted N^^ and jf^^ are reported for both the elec- 
tron and muon channels together with the Xc counting 
obtained with a single-gaussian fit where the Xc is con- 
sidered as a single peak. 

In order to verify the assumptions made, a system- 
atic study of the effect of releasing the different fixed 
parameters or varying them within a range around the 
assumed values was done, as well as a cross check of the 
Xc counting with the signal modelled as a single gaussian 
to describe both Xc states. The results of these studies 
are discussed in Sect. IVIIII 



VII. POLARISATION 

The experimental determination of polarisation can be 
used to probe assumptions on the impact of specific QCD 
processes and the infiuence of nuclear effects. The data 
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FIG. 3: /^"""/i invariant mass spectra in the muon channel for C (a), W (h), Ti (c) and full sample (d). The bin width is 
ISMeV/c^. 
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FIG. 4: e^e invariant mass spectra for C (a), W (h), Ti (c) and full sample (d). The bin width is 30 MeV/c 



Mat. 



H /i 



iV 



j/i, 



(MeV/c^) 



TV 



j/i/. 



(MeV/c^ 



C 80400 ± 300 35.6 ± 0.2 50030 ± 530 64.2 ± 0.8 

W 47750 ± 200 36.0 ± 0.3 23460 ± 480 66.1 ± 1.6 

Ti 4700 ±70 37.1 ±0.7 3530 ± 150 58.8 ± 2.8 

Tot 122900 ± 400 35.8 ± 0.1 77020 ± 700 64.3 ± 0.8 



TABLE VI: Total numbers of J/tp events per target material 
and for the full data set in the fi'^ ii~ and e'^ e~ channels. 



Mat. 


1-G fit 


2-G fit 


M^M" 




^Xc 


Nxal+Xa2 


^Xe2 


(GeV/c^) 



C 6280 ± 510 6390 ± 420 1.20 ± 0.26 0.030 

W 3120 ±560 2830 ±330 1.26 ±0.52 0.032 

Ti 390 ±110 390 ±110 0.63 ± 0.63 0.030 

Tot 9570 ±710 9630 ±550 1.19 ±0.24 0.031 



e+e- 


^Xc 


A^Xcl+Xc2 


'Vxc2 


(GeV/c^) 


C 3890 ± 480 


3600 ± 390 


0.79 ±0.31 


0.032 


W 2080 ± 370 


1870 ± 330 


0.71 ±0.48 


0.034 


Tot 5630 ±660 


5250 ± 500 


0.76 ±0.28 


0.033 



available for this analysis of Xc production does not allow 
a determination of the polarisation of the Xc states be- 
cause of large backgrounds. In the following we discuss 
angular distributions for the decay products of Xc and di- 
rectly produced J /"ip states with the goal of investigating 
the possible influences of polarisation on the acceptances 
and thus on the determination of the Xc rates. 



TABLE VII; Results of the fit of the Xc signal in the fj,^ j-i~ 
and e'^e~ channels. See text for the meaning of the different 
fit procedures. 
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FIG. 5: AM spectra in the muon channel for C (a), W (h), Ti (c) and full sample (d). 
background subtracted spectra the broken curves are the fitted Xci and Xc2 states. 
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FIG. 6: AM spectra in the electron channel for C (a), W (b) and full sample (c). The Ti spectrum is not shown due to the 
low statistics. The bin width is 10 MeV/c^. In the background subtracted spectra the broken curves are the fitted Xci and Xc2 
states. 



A. Xc polarisation 



isation. 



The full angular distribution of final state particles in 
the radiative decay 



XcJ -^ jJ/ip -^ 7 i~^l~ 



(6) 



can be found for pure Xc polarisation states | J, M) with 
J = 1,2 and \M\ = 0, ..., J in the appendix. The angu- 
lar distribution formulae are independent of the choice 
of a particular polarisation axis (e. g. Gottfried-Jackson, 
CoUins-Soper or other systems can be used). Possible 
coherent mixtures are not considered here because we 
assume that a study of the pure states will be sufficient 
to determine systematic acceptance effects due to polar- 



If one assumes no azimuthal dependence for the pro- 
duction process, the Xc decay depends on three angles 
which are chosen as follows: a polar decay angle, 9, defin- 
ing the direction of the J/ijj in the Xc rest system with 
respect to the polarisation direction; a polar angle, 9', 
defining the direction of the positive lepton in the J/ip 
rest system with respect to the J/ijj direction (in the Xc 
rest system); an azimuth angle, 4>' , which is the angle 
between the plane defined by the polarisation axis and 
the J/ip direction, and the decay plane of the J/tp. 

For a state \J,M) the angular distribution can be 
decomposed into terms with trigonometric expressions 



12 



and coefRcients K, 



J,M 



w 



J,M I 



b') = 



J2K''''Tfie,e',cp'). (7) 



The angular functions T/{9, 9', cj)') and the coefRcients 
Kj^ ' , expressed in terms of hehcity amphtudes, are re- 
ported in Table |X] in the appendix. With the additional 
assumption that for both Xc states only the leading mul- 
tipole, the electric dipole, contributes to the radiative 
decay, the coefficients K^ ' are uniquely defined (see ap- 
pendix for the numerical values). The assumption that 
higher order multipoles can be neglected is well justified 
by experimental results Q. The pure Xc polarisation 
states are thus unambiguously defined. 



B. J/tp polarisation 

1. J/ijj angular distributions 

In leptonic J/ip decays, the J/ip polarisation can be 
determined from the angular distribution of the leptons. 
After integrating over the azimuthal orientation of the 
decay plane of the J/ip (or assuming azimuthal symme- 
try) the distribution of the polar decay angle 6' can be 



parameterised as: 



1 dN 
Nd cos 9' 



a{X){l + Xcos^9'),a{X) 



1 



2(1 + A/3) 



(8) 



where 9' is the angle between the /''" and the quantisation 
axis. The form of the distribution ([8|) is independent 
of the chosen quantisation axis (in general however the 
value of A is dependent on this choice) . 



C. Method for the evaluation of systematic 
uncertainties due to polarisation 

In this section we explain the method used to estimate 
the systematic uncertainties for R^^ and R12 arising from 
possible polarisations of the Xc and directly produced 
J/tp states. The only experimental constraint which can 
be used for these estimations is the measured Xobs (sec- 
tion |VllB2|. 



1. Principle of the method 

The efficiencies entering in the formulae for R-^^ and 
R12 (ISl SI) depend in general on the polarisation of the 
Xc and the directly produced J/tp states. The efficien- 
cies will be evaluated for the Xc pure polarisation states 
described in section IVII Al which will then be used to 
limit the ranges of possible R^^ and R12 values which will 
in turn be used to determine the uncertainties of these 
values in section IVIIIBI 



As in the J/ip analysis |44| . 
we evaluate the polarisation states in the CoUins-Soper 
frame. Despite this specific choice and the restriction to 
pure polarisation states, we assume to get an estimate of 
uncertainties induced by polarisation. 

The formulae for R^^ and R12 require the detection 
efficency for direct J/%p production e* L which depends 
on the polarisation parameter Xdir- Since the observed 
polarisation, Xobs, also includes the effect of possible Xc 
polarisation, Xdir has to be disentangled from it using the 
values Al and A2 obtained for the assumed polarisation 
states of Xci and Xc2 respectively. This is done with 
an iterative procedure in which the yet-to-be-determined 
values of R^^ and R12 are used as inputs. 



2. J /ip polarisation measurement 

A measurement of J /ip polarisation by the HERA-B 
collaboration is reported in [4J]. However the J/tp sam- 
ple used for this study includes not only directly pro- 
duced J/f/) but also J l%p from Xc and it is not possible 
to distinguish between the two contributions. Therefore 
the A- value derived from the observed distribution, Xobs, 
has to be considered as the average J /ip polarisation pa- 
rameter, independent of the origin. 

The polarisation parameters have been determined us- 
ing as quantisation axis, the bisector of the angle be- 
tween pb and pt, where pf,, p^ are the momenta of the 
beam proton and the target nucleon, respectively, in the 
J/ip centre-of-mass system ('Collins-Soper frame'). The 
experimental value of the polarisation parameter, aver- 
aged over the muon and the electron decay channels and 
on the target materials, and assuming no dependence on 
Pj! and xj in the HERA-B acceptance, is |44l |: 



Xobs = -0.35 ±0.04. 



(9) 



2. Determination of Xdir, X\, X2 

Starting from Eq. ([S]) the observed polar decay angle 
distribution can be decomposed into contributions from 
directly produced J /ip and J/V' from Xci and Xc2 events: 

^ Y, ha^{l + X^cos'9'), (10) 

with Oj = a{Xi). The fractions fi{i = dir,\,2) of the 
different types of J/f/; are determined by i?^^ and Ri2- 
With EJ. - 1, /i + /2 = Rx. and /1//2 = R12 one 
obtains 



a, 



,hs(l + AobsCOS^^') 



Jdi' 



l-i?v 



/l 



RxcRi-2 

1 + R12' 



/2 



i?v 



I + R12 



(11) 



Since there is no direct measurement of Ai and A2, 
the angular distributions corresponding to the different 
pure polarisation states | J, M) of xa and Xc2 described 
in Sect. IVII Al are used to determine the sub-ranges al- 
lowed for Al and A2, out of the full [—1.0, 1.0] interval. 
For this purpose, the unpolarised Xc angular distribution 
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is re- weighted with the corresponding function ([7]). The 
resuhing cos 6*' distribution is then fitted with the func- 
tion ([8]) which then yields Ai or A2 corresponding to the 
tested pure Xc polarisation state. 

Solving (|10[) for Xdir as a function of Xo^s for given 
values of Ai, A2, Rx^ and R12 yields: 



A^' 



^diri^obs \ ^^17^2, RxaJ ^12) 



aobsXobs — O'l /l Al — 02 /2 A2 



flobs — a-i fi — a2 /2 

(12) 
In this equation R^^ and R12 enter via the fractions fi. 
On the other hand, as both depend also on Xdir, an iter- 
ative procedure is applied starting with Xdir = 0. 



3. Polarisation dependence of the efficiencies 

For each tested pure Xc polarisation state with the cor- 
responding set of values Xdir, Ai, A2, new efficiencies Sj 
and £jM are determined. 

a. J/ijj efficiencies: Assuming no dependence of A 

on Xp and Pj. (approximately valid within the un- 
certainty of our measurement |44|), we can write ej/^{X) 
as: 



ej/^(^) 



jSrJ/i> 
^ ^reco 



N, 



j/i, 



__ J A{cose',P) ■ M{P)- (1 + A cos2 61') -dcose'dP 
" / M(P) • (1 + A cos2 9')-d cos e'dP ' 

(13) 
where 9' is the polar angle in the polarisation frame, 
P is shorthand for all the other phase space variables, 
A{cos9',P) is the acceptance at the kinematical point 
{cos9',P), M{P) is the squared matrix element in P 
and A is the polarisation parameter. After calculating 
the integrals we find: 



£j/4,W = ej/^(^ = 0) 
where (cos^ 9') is given by: 
(cos^ 9') = 



2 a'\ 



1 -I- A • (cos 
(l + A/3) 



(14) 



2 _ J A{cos9',P) ■ M{P) -008^9' -d cos 9' dP 



J A{cos 9', P)- M{P)-d cos 9' dP 

All J/^ efficiencies, both for direct J/ip and for J/ip from 
the two Xc states, are calculated using Eq. [M) 

b. Photon efficiencies: To determine the effect of 
polarisation on the photon efficiencies we start with the 
formula: 






'7 ArXcJ.Af ' 



(15) 



where J, M denotes the polarisation state considered, 
N^7!i'.^' is the number of J/V' coming from Xcj and 



N^ 



'.//^ 



XcJ,M 
J/i, 



is obtained from a fit of the l~^l mass distri- 



bution, where each event enters with a weight: 

1 + Xjm ■ cos^ 9' 



'w{cos9',X,im) = 



(1 + Aj,m/3) 



(16) 



is the number of observed XcJ- The value of 



The value of Nx^ , m i^ obtained from a fit of the AM 
distribution of true XcJ (using MC generator informa- 
tion to select the correct J/ "07 combination), where the 
weight for each entry in the histogram corresponds to 
a certain pure polarisation state of XcJ, calculated by 
Eq.H 



VIII. SYSTEMATIC UNCERTAINTIES 

A. Uncertainties from reconstruction, calibration, 
simulation and background subtraction 

With the exception of J/ijj counting, all of the system- 
atic uncertainties in the measurement of Rx^ are com- 
mon to the e'^e~ and ijl^ ii~ channels. The J/^ counting 
systematic uncertainties are estimated to be 2% in the 
electron channel and 0.25% in the muon channel. 
The remaining systematic uncertainty estimates include: 

Xc counting: 

- photon selection (7%); 

- variation of I'^l^ mass window (2%); 

- Xc counting procedure (4%) including: 

• variation of the fixed parameters of the 
double-gaussian fit: AM^^^, cr^^i, AM(xc2) — 
AM(xci) and |^; 

• fit with free AM^^^ and/or a^^^ ; 

• change of binning of the AM spectrum. 

- Extensive tests were performed on the background 
determination with the mixed event procedure: 

• variation of corrections corresponding to com- 
binations of J/ip^s with photons from Xc de- 
cays in ME which do not occur in SE (+3%); 

• relaxing the requirement of the same neutral 
cluster muhiplicity in ME and SE (±2%); 

• variation of the cut on the neutral cluster di- 
rection in ME with respect to SE (±3%); 

• allowing for an additional, polynomial term 
in the background to improve the fit around 
the Xc signal yields an asymmetric uncertainty 

(+4%). 

The total contribution from the background de- 
scription to the systematic uncertainty in the Xc 
counting is estimated to range between —4% and 

+6%. 
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Efficiency evaluation 

- the use of different kinematic distributions for the 
generation of the J/i/" (see Ref. [28] and [STi]), af- 
fecting both "-i^ 
effect on R^ 



and £^, introduces a systematic 

'■■j/'p 
of 4%; 



- tests on the photon efficiency simulation were per- 
formed including the comparison between real data 
and Monte Carlo of the 7 conversion yield and of 
the detection efficiency of photons from electron 
Bremsstrahlung. The overall systematic uncer- 
tainty on e^ determination and correction is found 
to be 6.5%. 

The overall systematic uncertainty on R^^, evaluated 
as the quadratic sum of the above terms, is therefore 
1^2 "^ for both J/ijj decay channels. 

The systematic uncertainty on R12 is completely dom- 
inated by the accuracy of the ECAL energy calibration 
which affects AM(xci) which in turn affects the ratio 
J^"^ . A fine tuning of the ECAL calibration as a function 

of the photon energy was performed using the tt'^ — > 77 
signal. An absolute calibration accuracy of ^ 2% was 
obtained and on AM^^^ of ^ 8 MeV/c^. By scanning 
AM{xci) in such range around the nominal value [5|, a 
variation of j^^^ (and thus R12) of 35% is obtained. No 
effect on R12 is observed by changing the other fixed fit 
parameters. 



B. Polarisation effects 

Since the direct J/ip and Xc polarisations cannot be de- 
termined separately from our data, we estimate instead 
systematic uncertainties on the reference values reported 
in Table Em] (and denoted by R^,''/ and R[l^ in the fol- 
lowing) which were obtained with the assumption of zero 
polarisation. The results of this study are expressed as 
overall shifts of the values of R^^ and R12 due to the 
average polarisation of directly produced J/^p with un- 
certainties obtained from the maximum variation of Xc 
polarisations allowed by the measurement: 



Rti^—RZ. 



R12-R 



-9.5%tv:" 



(17) 



where the following ingredients are used: 

• The central values are obtained from the average 
measured value for Aobs and with the assumption of 
no polarisation of Xci and Xc2 (Aofcs = —0.35, Ai = 
and A2 = 0, yielding Xdir = -0.424). Therefore, if 
the observed J/^ polarisation were due exclusively 
to direct J/ip polarisation, the measured R^^ would 
be shifted up by 9.5%, while obviously no effect on 
R12 is produced. 



The variation bands in Eq. [T7| are obtained by tak- 
ing the extreme positive and negative variations of 
the central values defined above, of all combina- 
tions of Aofcs (varied in a 95% c.l. range around the 
measured value, see Eq. ^ with Ai and A2 (cor- 
responding to the different pure helicity states Ml 
and M2): 



upper value: Aobs = 
0.18 for i?x^; Ai = 

lower value: Xots 
A2 = -0.18 for i?x 
for i?i2. 



-0.44; Ai = -0.24, A2 = 
-0.24, A2 = -0.18 for E12; 

= -0.26; Ai = 0.22 and 
; Ai = 0.22 and A2 = 0.18 



• Different polarisation values give overlapping 
ranges of possible R^^ and i?i2 values. Any value 
in each range is equally probable. Thus, even if 
the error on Xots was Gaussian distributed, the er- 
rors of R^^ and R12 would not be Gaussian dis- 
tributed. To take into account that the polarisa- 
tion parameter Xots was determined as an average 
over the whole accepted phase space and over dif- 
ferent materials, Xobs was varied in a ±2(7 range 
with equal weights. Selecting the maximum devia- 
tions the measured values i?^^and R12 would have 
to be scaled: 



^xc = /i?x<= ■ ^x/ 



R12 — lui-i ■ Rr2 ' 



(18) 



with JR^^ e [1.02, 1.21] and Jr., e [0.89, 1.16], 
where the uncertainties due to polarization are fully 
contained in the ranges given. 

Note that the correlation between the values of R^^ 
and i?i2 are ignored in Eqs. [TTlandfTSl 



IX. RESULTS 



R\ 



The measured values for R^^ are computed from Eq. [31 
assuming zero J/ip polarisation and are reported in 
Tab. IVIIIi separately for muon and electron channels and 
combined sample. When averaged over decay channels 
and target materials, a value of 



Ry 



0.188 ±0.013^ 



h0.024 
-0.022sj^s 



(19) 



is obtained. The quoted uncertainties include all system- 
atic contributions (except the polarisation contribution 
which is given in Eq. [T7]as a variation band at 95% c.l). 
The following observations can be made: 

• The results obtained in the two lepton channels are 
compatible within la in both G and W samples. No 
measurement for the Ti in the electron channel is 
possible due to the low statistics; 
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• The values of R^^ obtained separately in the three 
target samples are consistent with each other; 

• The present result is lower than most values pub- 
lished in the literature in pN interactions (see 
Tab, mil and Fig. [T]). Despite the fact that the vari- 
ous available measurements are taken at widely dif- 
fering centre of mass energies, they are for the most 
part compatible within r^ l.Scr, except for E705 
(2.3cr) andlSR (3.3cr). 

The present measurement is lower than the pre- 
vious HERA-B resuh [li| by about 2a. The two 
analyses are quite similar, although more extensive 
systematic checks have been performed in connec- 
tion with the present one. These checks did not 
uncover any error in the previous analysis and we 
thus believe that the differences are largely sta- 
tistical. The average of the two HERA-B results, 
R^^ = 0.198+0026' differs by less than la from the 
result of Eq. [HI 



B. 7?i2 

The measured values of R12 are evaluated using Eq.[31 
assuming no polarisation for either the directly produced 
J/tjj^s or the Xc's, and are summarised in Table FVllIl 

As above, no dependence on target material is ob- 
served. The results from the electron channel are con- 
sistently lower than the muon results, but nonetheless in 
agreement to within la of the statistical uncertainties. 

The final result averaged over decay channel and target 
material is: 



i?i2 



1.02 ±0.17,, ±0.36 



sys 



(20) 



where the systematic uncertainty does not include the 
polarisation contribution which is given in Eq. [T7| as a 
variation band at 95% c.l. The J/V' yields from Xci and 
Xc2 are therefore found to be equal, although with large 
uncertainties. 



C. Dependence on kinematic variables 



A study of the dependence of R^^ on the kinematic 



J/^ 



J/^ 



variables Xp and pj, in the ranges covered by HERA- 

B {xy^'' e [-0.35, 0.15], p^/'^ < 5 GeV/c) was performed 

by applying the described procedure in five xj and 

three pj intervals respectively. The resulting distribu- 
tions, for both channels combined, are shown in Fig.[7^)- 
b). The data is compatible with a flat dependence of R^^ 
on both kinematic variables, although more complex de- 
pendences cannot be ruled out. 



D. A-dependence 

The atomic mass number (A) dependence of inclusive 
cross sections is often parameterised as a power law: 



CpA = CTpNA" 



(21) 



where apj^ is the inclusive production cross section in 
collisions of protons with a nuclear target of atomic mass 
number A, apN is the average cross section in collisions of 
protons with a single nucleon and a characterises the A 
dependence of the cross section. The difference between 
a for J/ip production and that for Xc production can be 
computed from the measured values of R^^ for C and W 
targets given in Tab. IVlITl from the following formula: 



Aa 



I. 



ay 



'J/'4' 



log^ 



-ftv 



W 



c 



(22) 



where Aw — 184 and Ac = 12 are the tungsten and 
carbon atomic mass numbers. The results, plotted as 



a function of 



Jl-<i> 



and p-p 



are shown in Figs. [Th), 



Averaged over the visible 



Jl^ 



and p'^ range. 



d). 

Aa = 0.05 ± 0.04. The predictions of the various pro- 
duction models for Aa are all within the uncertainties of 
the measurement 1231. 



E. 



Xc 



cross sections and ratio 



From Eq. [5] we obtain the values for the cross section 
ratio '^^^^ \ under the assumption of zero polarisation for 
both J/^ and Xc- The results are reported in Tab. IVIlH 
The target material averaged result is: 



o-(Xci) 

0'(Xc2) 



0.57 ±0.23 



(23) 



where the uncertainty includes the systematic contribu- 
tions (except polarisation - see above). The Xc produc- 
tion cross sections, defined as: 



cr(Xc 



a{J/^)R^ 



Br{xc^ -^ J/ipl) ' 



1,2 



(24) 



are calculated using the estimate of the total J/tp 
cross section at -^/s = 41.6 GeV, a{J/^p) = (502 ± 
44) nb/nucleon reported in [29| and assuming that R^^ 

is independent of Xp over the full xj and prp range. 
The following target material averaged values are ob- 
tained: 



""(Xci) — 133 ± 35 nb/nucleon; 
c(Xc2) = 231 ± 61 nb/nucleon. 



(25) 



leading to a total Xc production cross section cr(xc) = 
364 ± 74 nb/nucleon. Fig. [T] shows all available measure- 
ments of the Xci and Xc2 production cross sections and 
their ratio in proton-nucleus interactions at fixed-target 
energies. 



Mat. 



C 
Ti 
W 
Tot 



MM 



combined 



Ry 



0.n4:±0.029sttQ°oll, 



0.190 ±0.018,tt""^* 

ys ^-^^-^ sys 

+0.025 



0.185 ±0.015at 1^022, 



0.197 ±0.079, 



s* -0.023 sj/s 



0.197 ±0.079, 



+0.025 
*~0.023sas 



0.202 ± 0.055st taoM,^^ ^-^^^ ^ ^■^'^'^st. 0,022 sys 



0.191 ±0.034, 



0.180 ±0.025,(1^2? 



0.190 ±0.015, 



0.194 ±0.029,tl:o:o23, 
0.188 ±0.013,t+n^^ 



Ri 



C 
Ti 
W 
Tot 



0.82 ±0.32, 

0.73 ± 0.49, 
0.79 ± 0.27, 



1.23 ±0.27,4 1.06±0.21,t ±0.37,„, 

0.67±0.67,t 0.67 ±0.67,t ±0.23,,;, 

1.27±0.53,t 0.98 ± 0.36,t ± 0.34,„, 

1.17±0.22,t 1.02±0.17,t ±0.36,„, 



°"(Xcl) 
a"(Xc2) 



c 

Ti 

W 
Tot 



0.47±0.19,t 

0.41±0.28,t 
0.45±0.16,t 



0.70±0.16,t 0.60 ±0.12,t± 0.21, 

0.38±0.38,t 0.38 ±0.38,t ±0.13, 

0.72±0.30,t 0.56 ±0.21,t ±0.20, 

0.66±0.13,t 0.57 ±0.10,t ±0.20, 
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TABLE VIII: Measured values of R^^, R12 ant 
full data sample. 



g"(Xcl) 
ct(Xc2) 



/i"*"^ and combined sample for the different materials and the 



u 0.4 

?^ 

0*0.35 


— 
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= — 1 
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0.15 
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0.1 
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0.05 


1 








-0.35 



^ 



-0.085 
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FIG. 7: Dependence of R^^ on xj (a) and Pj! (b) (circles: C; triangles: W). Dependence of Aa = a^^ — aj/^ on xj (c) 
and p'j! (d). Only statistical errors are shown. 



X. CONCLUSION 

We have presented a new measurement of the fraction 
of all J/'0 mesons produced through Xc decay (i?^^ ) , per- 
formed with the HERA-B detector in pC, pTi and pW 
interactions at 920 GeV/c (^s = 41.6 GeV). The Xc 
mesons were detected in the J/V'7 decay mode, and the 
J/ij) in both /i^/i^ and e+e^ decay modes. The de- 



J/^ 



and extended from 



tector acceptance was flat in prji 

x-^J"^ = -0.35 to Xp^"^ = 0.15. 

The measurement is based on a total sample of ~ 15000 
Xc, the largest ever observed in pA collisions. Apart from 
lepton identification requirements, the analysis is identi- 
cal for the two channels. The separate results for the two 
channels are found to be in agreement with each other in 
all respects. 



The measured value R^^ ^ 0.188 ± 0.013,tt[J:[]^^^ ^ is 
~ 2a lower than the previously published result from 
HERA-B. Our new value is also lower than, but not in- 
compatible with, most of the previously published val- 
ues obtained from pN interactions, independent of the 
centre of mass energies and the kinematic ranges of the 
measurements. The present result supports the NRQCD 
calculations [23]. When taken together with the already 
published result of HERA-B on ijj' production [28], the 
fraction of all J/4> mesons coming from decays of higher 
mass charmonium states is found to be ~ 27%. 

By separately counting the contribution of Xci and Xc2 , 
wc obtain a ratio of the two states R12 = ^Xci/^xc2 ~ 
1.02 ± 0.40 and a cross section ratio 4^ = 0.57 ± 
0.23. The Xci and Xc2 cross sections are measured to 
be cr(xci) = 133 ± 35 n6/nucleon and (t{xc2) = 231 ± 
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61 n&/nucleon in the full xj range. 

No significant departure from a flat dependence of R^^ 

on the kinematic variables xJ and p-J is found within 
the limited accuracy of our measurement. No significant 
difference in the A-dependence of Xc and J /ip production 
is found within the limits of the available statistics. 

For the first time, an evaluation of the effect of po- 
larisation of J /ip and Xc on the measured values of R^^^ 
and i?i2 was performed. The behaviour of R^^^ and R12 
as a function of the polarisation, expressed by the A pa- 
rameter, was studied with the conclusion that R^^ and 
R12 are uncertain with factors in the ranges [1.02,1.21] 
and [0.89,1.16], respectively, ignoring correlations be- 
tween the two. 

No mention of the influence of polarisation on the mea- 
surement of Ry^^ can be found in any of the previous 
measurements. Nonetheless, we suspect that all mea- 
surements are subject to similar uncertainties to greater 
or lesser extents, depending on the geometry of the ap- 
paratus used. 



with the density matrix for the J/^ helicity (jIX 



(A2) 



K = ±l 



Using the notation of [43| the angular distribution can 
be decomposed into terms with trigonometric expressions 
Tf{e, 9', (j)') and coefficients K--^^\A^^^): 



W'-''{e, 6', 0') ^Y^K^^'^^A^i^) T/{e, e\ cb'). (A3) 

The Kf-^^{A'^ and T/{e, 6', (j)') are reported for J = 
1,2 in Table 1x1 The normalisations are for the angular 
distributions 
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APPENDIX A: Xc ANGULAR DISTRIBUTIONS 



w'^^^{9, e\ (j)') d cos ed(i>d cos e'< 



647r^ 
9 



for J = 1 



^ for J = 2 

15 

(A4) 



The helicity amplitudes A'/^ can be expanded in multi- 
pole amplitudes (El, M2, E3), see for example ,43;]. With 
the restriction to electric dipole transitions the helicity 
amplitudes become: 



J= 1 : 



An 



Ai 



(A5) 



J — 2: Ao — \/-rT:, Ai 



A^ 



The angular decay distribution of a pure polarisation 
state \J,M) is given as an expansion into helicity ampli- 
tudes ^/ I bv I43I : 



W 



J,M I 



h 0', cj,') = 



(Al) 



+J 



I/, y'^ — J /^— ±1 



The TablelXlreports also the coefficients K^ ' calculated 
with these values for the helicity amplitudes. Hence, with 
the restriction to the lowest multipole, the angular distri- 
butions of a Xc decay for a given polarisation state | J, M) 
is fully determined (obviously, the relative contributions 
of different polarisation states are not fixed) . 
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